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A METHOD FOR PREDICTING SHOCK SHAPES AND PRESSURE 


DISTRIBUTIONS FOR A WIDE VARIETY OF BLUNT 
BODIES AT ZERO ANGLE OF ATTACK 
By George E. Kaattari 
Ames Research Center 


SUMMARY 


A method is presented for determining shock envelopes and pressure dis- 
tributions for a variety of blunt bodies at zero angle of attack. Correlation 
functions obtained from exact solutions are used to relate the shock standoff 
distance, at the stagnation and sonic points, to the body geometry. These 
functions were obtained for a perfect gas but may be applied for real gases in 
equilibrium flows. The method is restricted to cases where the bow shock is 
detached from the body and the flow over the forward face is subsonic. Results 
given by the method are shown to be in good agreement with experimental values. 


INTRODUCTION 


A problem currently receiving wide attention is that of predicting the 
shock envelope and pressure distribution over the forward face of axisymmetric 
blunt bodies at atmospheric entry. 

Exact solutions have been presented for spheres, ellipsoids, and parab- 
oloids by Van Dyke (ref. 1). In this work and that of subsequent investiga- 
tors most exact solutions are for the indirect problem; a shock shape is pre- 
scribed and the resultant body shape is sought. Many trial solutions are 
usually necessary to find a body shape that approximates the one for which a 
solution is wanted. 

An approximate "direct” method to define the shock trace at zero angle 
of attack is included in the analysis of reference 2. This method is 
restricted to spherically blunted bodies with sharp corners. Extensions of 
this method which permit the calculation of the entire forward shock envelope 
and the pressure distribution for these bodies have been published (refs. 3 
and 4 ). The exact and approximate solutions of these references are applicable 
to a relatively narrow choice in body shape. 

This investigation was undertaken to develop methods for predicting the 
shock envelopes and pressure distributions over the forward face of a wide 
variety of body shapes at zero angle of attack. Particular attention is 
given to ellipsoids (ranging from a sphere to a flat disk), spherical or 



flat -faced bodies -with an arbitrary degree of corner rounding, and spherically 
blunted, large -apex-angle cones. The method developed is restricted to those 
cases in which the bow shock is detached from the body and the flow over the 
forward face of the body is subsonic. 


NOTATION 


A coefficient in transformed pressure distribution function 

a semiminor axis of ellipsoid 

body bluntness parameter, y 2 = 2R b X] 3 - B^x^ (ref. 1) 

B s shock bluntness parameter, y 2 = 2R s x s - B g x| (ref. 1) 

b semimajor axis of ellipsoid 

C coefficient in transformed pressure distribution function 

G shock correlation function (see eqs. (1)) 

M free -stream Mach number 

m* slope of transformed pressure distribution function at stagnation 

S point 

m£ slope of transformed pressure distribution function at sonic point 

p pressure on body 

p pressure at stagnation point on body 

s "C 

p^ pressure at sonic point on body 

R^ body surface radius on axis of symmetry 

R^* body surface radius at sonic point of an ellipsoid 

^b§ "body surface radius at point on surface inclined $ degrees with 

respect to a plane normal to free -stream direction 

R_ shock radius on axis of symmetry 

o 

r c corner radius of body 

r m semidiameter of body 

s distance along body surface from stagnation point 
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distance along body surface from stagnation point to sonic point on 
body 

independent variable for transformed pressure distribution function 
tangent point of transformed pressure distribution functions 
free-stream velocity 
velocity at sonic point on body 

streamwise distance from apex of body to sonic point on body 

streamwise distance from apex of shock to point on shock at distance 
y^ from axis of symmetry 

normal distance from axis of symmetry 

normal distance from axis of symmetry to sonic point on body 
coordinate axis of transformed pressure distribution function 
dependent variable in transformed pressure distribution function 
specific heat ratio 

shock standoff distance from body in free-stream direction 

shock standoff distance from body on axis of symmetry in free-stream 
direction 

shock standoff distance from sonic point on body in free-stream 
direction 

inclination on forward body surface at the tangent point with corner 
radius^ r c 

normal distance from body surface (sketch (g)) 

shock surface inclination at distance, y, from axis of symmetry with 
respect to plane normal to free-stream direction 

shock surface inclination at a point opposite the sonic point on the 
body with respect to a plane normal to the free-stream direction 

shock surface inclination, at a point opposite the sonic corner of a 
flat disk, with respect to a plane normal to the free-stream 
direction 

pressure correlation function 
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p 1 density of free stream 

p 2 density of stream behind normal shock 

p s _k density of gas at stagnation point 

o shock layer thickness normal to body (sketch (g)) 

0 body surface inclination with respect to plane normal to free -stream 

direction 

body surface inclination at sonic point with respect to plane normal 
to free -stream direction 

0^ body surface inclination of a sphere at the sonic point with respect 

to a plane normal to the free-stream direction 


ANALYSIS 


The analysis will be developed in three parts. First, general relation- 
ships and assumptions involving body and shock geometry are introduced. Next, 
application of the shock-body relationships to a method for calculating the 
shock envelopes of a variety of blunt bodies is demonstrated. The calculation 
involved is simplified by the use of a nomograph. Finally, a general method 
for estimating the pressure distribution for the bodies considered is 
presented. 


General Shock -Body Relationships 

Shock shape . - Reference 1 defined shock traces as conic sections of 
suitable bluntness B s (see Notation) to obtain the shock solution for a 
range of conic section blunt bodies. In the present investigation, it was 
assumed that conic sections can closely approximate the shock trace for a 
larger variety of blunt bodies than were considered in reference 1. 

Stagnation-point cor re latio n. - A relationship involving shock standoff 
distance, shock shape, and normal-shock density ratio was pointed out in 
reference 5 and used in reference 2. A simplified form of this correlation, 
as used in the present investigation, is 



(la) 


(lb) 
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G is a function of the free -stream gas and Mach number and was evaluated by- 
substituting values for A 0 /R s and A 0 /R- b from reference 1 into equation (la). 
Figure 1 shows G as a function of the reciprocal normal-shock density ratio 
p 1 /p 2 (rather than Mach number) for two values of the gas specific heat ratio, 

7- 


Son ic point shock standoff correlation . - The theoretical results of 
reference 1 and numerous experimental results for a variety of body shapes 
indicate that the ratio, A^/y^, is primarily a function of the normal-shock 
density ratio and secondarily of y pertinent to the free -stream gas in ques- 
tion. Sketch (a) shows A^, the shock thickness measured in the free -stream 
direction from the sonic point on the body at y*, the vertical distance from 
the axis of symmetry. 

A correlation of both theoretical and experimental values of A^/y*. with 
the normal -shock density ratio, p^/pg* is shown in figure 2. The dashed 
extension of the y = 1.4 line is physically unreal, but it is a useful inter- 
polation guide for A */y* at values of y between 1.0 and 1.4. 



Figure 1.- The G function. 


Sonic point inclination . - A 
solution for the infinite shock- 
density ratio was used to determine 
the inclination of the sonic point 
on conic-section bodies. The 
results were then applied to finite 
shock-density ratios and to bodies 
of other shapes. 

The angle, $*, associated with 
the sonic point is the angle between 
a normal to the free -stream direc- 
tion and a line tangent to the body 
at the sonic point (sketch (a)). 

This angle was determined with the 
Busemann solution for an infinite 
shock -density ratio as a function 
of the body bluntness parameter. 



Sketch (a) 
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B b (see Notation). Details of the solution are presented in appendix A. 
essential result is 

[(B b - 1 ) s in 2 $*+ l] 2 


2 (B b - 1 ) sin 4 >* 


The 


f 1 . 

J Bt, - 1 sin 3 >* 

sin 

[•K- 1 

«/( B b - 1) sin 2 $*+ 1 

- l)sin 2 + 1 

J 


i.o 


Body shape 


Symbol References 


Spherical (round cornered) O 4 

.9|— Sphericol (sharp cornered) k Unpublished 

Disk (round cornered) 

Disk (sharp cornered) 

Cone ( € =70°) 

,8|— Sphere 

Parabojoid 

] — Experiment - open symbols 

Theory - solid symbols 


D 6 

□ 4 

A Unpublished 



Figure 2.- Correlation curves of shock-standoff 
distance at sonic point, A 


= (l - sin 2 $„) - (2) 

Pst 

In the present application of 
equation (2), the gas was assumed 
perfect. Since the theory requires 
an infinite shock-density ratio, the 
value P*/P s t = 0.607, appropriate 
to a gas with 7 = 1.0, was used. 

The results of equation ( 2 ) for 
sin are presented in figure 3 
normalized with respect to the value 
for a sphere (sin 0 ^) . This nor- 
malized value is plotted as a func- 
tion of the more tangible variable, 
a/b = B b 1/2 , rather than as a 
function of B^. 

Although the absolute value for 
sin <£>* given by equation (2) may 
not be correct for a finite density 
ratio, it is assumed that the ratio 
sin <k*/si n will not differ 

greatly from that given by exact 
theory. 1 Thus, if the value of 
sin is known for a sphere in a 

given free stream, the value of 
sin ^ for an ellipsoid can be 
determined from the curve of fig- 
ure 3 * Values of for a sphere 

at various free -stream conditions 
(ref. 1) are presented in figure k 
for convenience in such 
calculations. 

If a blunt body does not have 
a conic section, the parameter B b 
is not appropriate. However, the 

^■The validity of figure 3 for 
finite shock-density ratios is sub- 
stantiated to some degree by the 
points from reference 1 for ellip- 
soids of moderate values of a/b. 
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sonic -point location on other types 
of blunt bodies may be estimated on 
the basis of the foregoing analysis. 
To this end, it is assumed that a 
shock in the form of a conic section 
is appropriate to all classes of 
blunt bodies considered here. For 
example, a blunt body with a rounded 
corner of radius r c is illustrated 
in sketch (b) . When the flow has an 
infinite shock-density ratio (p 2 /p 1 )^ 
the shock approaches coincidence 
with the body in the vicinity of the 
sonic point. If the sonic point is 
on the rounded corner of the body, 
the above Busemann solution applies 
if the body corner radius r c is 
equated to the radius of a 

conic-section body at the same sonic 
point ordinate, y^. 


Figure 3.- Sonic angle, on ellipsoids. 


R b* 


sin $*[ (Bb - l)sin 2 <1 + 1] 



Figure 4 .- Sonic angle for a sphere 
(Bj, = 1, a/b = 1). 




sin 0^) 


or 


r 

r 


1 

1 + (B^ - 1) sin 3 ^ 


(3) 


Equation (3) gives the r c /r m 
value for an equivalent body that 
supports the same elliptical shock 
and the same sonic angle as a body 
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with the bluntness parameter, B^. The relationship between B b and 0^. was 
previously derived and is given by equation ( 2 ). The value of sin ^ nor- 
malized with respect to the value for a sphere is plotted in figure 5 as a 
function of r c / r m as determined from equation (3). The general validity of 
this function at finite shock-density ratios will be discussed in a later 
section. 



Figure 5-- Sonic angle, 

bodies. 


on round -cornered 



Sketch (c) 


Figure 5 is intended to apply 
to any blunted body with a rounded 
corner. In many cases, however, a 
qualification of the sonic angle is 
necessary, for example when the body 
is conically or spherically blunted. 
The conical or spherical portion of 
the body will be tangent to the cor- 
ner radius, r Q , at which point a 
surface inclination angle, €, is 
defined with respect to a plane nor- 
mal to the body axis of symmetry. 

The sonic angle is then taken as e 
if e is greater than 
determined from figure 5* 

Relationship between shock and 
body inclination. - Reference 2 shows 
that in a given free stream the 
shock angle, 0*., opposite the sonic 
point on the body is uniquely 
related to the sonic point angle, 

0 ^. The orientation of these angles 
is shown in sketch (c). The results 
of reference 2 are restricted to 
sharp -cornered bodies and are based 
on a gross mass 'flow continuity 
analysis. It is nevertheless 
assumed that such a unique angle 
correspondence does exist and holds 
for all bodies with a detached shock. 
If this assumption is valid, the 
functional relationship between 0 ^ 
and can be established from any 


set of shock solutions. The solutions of reference 1 for conic-section bodies 
were used for this purpose in the following manner: The value of y / R Q was 


found by solving equation ( 4 ) (appendix B) ; values of B b , B s , R s /R]3j and 
A q /R s were obtained from reference 1 and the value A^/y^ from figure 2: 


^*y* _ ^o 




■x-^s R S 


R> 


*b R s 


1 - 


!± - B b 


R s y* 

R b R s- 


_1_ 

Bo 


B r 


y 


* 


R c 


( 4 ) 
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With the value of y_^/R s determined, the following equations were then solved 
for 0 ^ and 0 ^: 


tan 9 * 



tan 0 ^ 


y* R s 

R s R b 



(5) 


( 6 ) 


The above procedure, giving the relationship between 0 ^ and was 
restricted to the limited range of body bluntness (B^ < 3 ) (as found in 
ref. 1) and, therefore, to a limited range of Additional procedures were 

necessary to determine the value of 0 ^ at smaller values of associated 

with bodies of larger bluntness. A simple solution for 0 ^ in the limiting 
case of a flat disk (B^ = °°, $* = 0°) was found (appendix C). The relation- 
ship between 0 ^ and 0 * over the range of 0 ^, not amenable to available 
theory was then approximated in the following manner. 

It was assumed that if the corner of a flat disk was rounded to a small 
radius, no important change in the shock shape would occur . 2 For the sonic 
point on the rounded corner the inclination could be significantly dif- 

ferent from zero. In the limit of a vanishingly small corner radius, no 
change occurs in the shock and, therefore, d0^. =0. It then follows that 
d 0 */d<i>^ = 0 at ^ = 0 °, suggesting that 0 ^ would develop in an even power 
series in 0 ^. Accordingly, a three -term power series in 0 ^ was used to 
approximate the variation in 0 ^ from its value at = 0 ° to the values of 
0 ^ (and d 0 ^/d^^.) in the range of 0 ^ determined from the results of 
reference 1 . 


2 The small variation in shock shape with changes in body shape is 
pointed out in reference 1 . 
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(a) as a function of 

Figure 6.- Inclination angles and 0^. 

Plots of 0^_ as a function of 0^ are presented in figure 6 (a) for 
various shock-density ratios and 7 values of 1.4 and 1.0. A useful corre- 
lation between the results for the values of 7 is presented in figure 6 (b). 
The ordinate represents the incremental value or change in Q (defined as 
0* - 0* o ) as a function of 0 ^. At a given shock-density ratio, this 
incremental change in 0^ is independent of 7 . 

These relationships between and do not completely agree with 
those of reference 2 , primarily because of the restrictions imposed on the 
shock shape in the earlier work. 

The manner in which the 6 ^ functions are applied in determining 

shock solutions for the classes of bodies considered here will be discussed 
in the following sections. 
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Shock Envelope Solutions 


Shock solutions will be considered for the following classes of bodies: 
ellipsoids ranging from a sphere to a flat disk, spherical or flat -faced 
bodies with an arbitrary degree of corner rounding, and spherically blunted, 
large -apex-angle cones. Equations relating the shock shape with the geometry 
of these bodies will be given. The use of a nomograph which simplifies the 
solution of these equations will be discussed. 


Shock equations for conic -section body.- A conic-section body of known 
bluntness, B^, and a shock bluntness, B s , at a distance, A Q , from the body are 
depicted in sketch (d). The geometrical relationships are: 



(7c) 


Equation (7c) is divided by equation (7a) and the result is rearranged to give 
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Non-conic -section body shock equations . - A section of a body whose 
geometry 3 involves all the elements of the non-con4c-section bodies to be 
considered is shown in sketch (e). The geometrical relationships are: 


y* = r m - r c (l - sin 



(9a) 


Xi 


: b = r m tan e + R b 


COS € - 1 
COS € 


. , 1 - sin e , 

+ r ^ [ - cos ® 

COS € * 


x s + A * = A o + x b = A o + r m tan e + R b 
1 - sin € 


+ r r 


- COS ® 

COS € * 


(9b) 

N 

C OS € - 1 

cos e 
(9c) 


Sketch (e) 

Equation (9c) is divided by equation (9a) and the result is rearranged to give 


Ad _ A 

y * y * 


r m tan e + 


c os e - 1 
cos e 


+ r, 


1 - sin e 
cos e 


- cos 


# 


x r 


r m - r cd - sln **) 


or 


Rg 

y* 


R> 


A* 

y* 


tan € + — — 
m 


L b f cos e - 1 


cos e 


m 


1 - sin e 
cos € 


- cos <£> 


1 (l - sin 0 ) 

T * 


A„ 


fs 

y* 

Aq 

R c 


( 10 ) 


Shock solution nomograph * - The equation of a shock, with the ordinate set 
equal to y^, is 


B s 5s. = ( h. - 


y* \y* tan 0. 


( 11 ) 


3 An Important additional restriction of the method applies to spherically 
blunted, large -angle cones: The cone angle (90° - e) must be sufficiently 

large to assure a detached shock even if the cone apex is sharp (R b = 0) . 
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Since the body geometry is given, can be determined from figure 3 or 5 and 

the corresponding value for Q can then be found from figure 6. Equa- 
tions (l), (8) (or (10)), and (11) can be solved simultaneously to yield the 
shock solution in terms of y^/R g , B , and Z^/Rg- This solution requires con- 
siderable computation. The computation can be minimized by use of a nomograph 
consisting of constant B s and 6 * curves on an R s /y* versus x s /y^ coordi- 
nate system. Such curves, calculated with equation (ll) and plotted in 
figure ar e used as follows. 



y * 


Figure 7«- Shock-solution nomograph. 

The numerator of the bracketed term in equation (8) or (10) contains 
elements of known value. If the value of A 0 /R s were also known, equa- 
tion (8) or (10) could then be plotted as a straight line on the coordinate 
system. This straight line would have the ordinate value equal to the 
bracketed term at x s /y^ = 0 and would have a slope equal to R g /A 0 . The 
shock solution would then be given by the intersection of the line with the 
nomograph curve having the appropriate known value of 0 # . This intersection 
point simultaneously locates the correct values for B g and R s /y*. The shock 
shape and its location with respect to the body are then completely defined. 




In most cases, Aq/R s I s not known so the procedure described above can- 
not be applied directly. However, a quick convergence to a solution is possi- 
ble through an iterative technique beginning with a gross initial assumption 
for the value of Aq/R s . Details of this iterative technique will be 
clarified in the Numerical Examples section. 


Pressure Distribution 


The method for estimating the pressure distribution on blunt bodies 
involves finding the simplest pressure distribution curve that conforms to 
"end" conditions specified at the stagnation and sonic points on the body. 
A suitable, easily determined curve was obtained through the use of the 
following coordinate system. 


Coordinate system . - In the natural coordinate system of pressure and 
location points on the body, pressure can be expressed as some power series 
in terms of distance from the point of maximum (stagnation) pressure on the 
body. Experiments indicate that for a flat disk the pressure gradient is 
infinite at the sonic point (sketch (f)). For such a case, an infinite power 


/ 



series in terms of distance from the sonic 
point is required to define the pressure 
distribution. In the present analysis, 
only a limited number of derivatives could 
be specified for the pressure distribution 
so that a term-by-term evaluation of an 
unlimited power series was not possible. 
However, if a pressure distribution curve 
without inflections could be represented 
in a coordinate system by a curve whose 
maximum gradient did not exceed a rela- 
tively small value, it would then be pos- 
sible to approximate such a curve with a 
simple expression of few terms. The coor- 
dinate system in sketch (f) was con- 
structed as follows: The origin of the 

coordinates was at s/s^. = 0 and p/p ^ = 1- 
The axis of the independent variable, 
s’/s^*, was directed from this origin to 


the point p/p 


st 


p */ p st at s / s * 
ate 


= 1 . 


sketch (f) The transformed coordinates are indicated 

by dashed lines. The transformed vari- 
ables are s '/ s ± and z. The absolute value of the slope of the transformed 
pressure distribution curve, z = f(s T /s^), cannot exceed the finite value, 
m, 


* 


= l/m^, as indicated. At s T /s^ = 0 and 1, z = 0 


The variation of z with s */ s * was approximated by two curves having 
appropriate slopes at s' /s£ = 0 and S V S * = respectively, and which were 
matched to a common tangent. The evaluation of these slopes, the resulting 
transformed pressure distribution functions, and the method of matching these 
functions to a common tangent are described in the following paragraphs. 
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Stagnation poin t.- The relationship between pressure and velocity for a 
perfect gas is 


7 


p _ 

fi -7 - 1 fxY 1 

Pst 

> 

> 

i—l 

+ 

: 


which, for small velocities, V, near the stagnation point, can be written as 


P 


P 


st 


£ 1 


7 

7+1 



( 12 ) 


The velocity variation with distance, s, is linear near the stagnation point 
so that equation ( 12 ) may be written, in natural coordinates, as 


JE_ = i 

Pst 



(13) 


The bracketed term in equation ( 13 ) is discussed in detail in appendix D. The 
transformed function, z, in equation (l4) is obtained from the Taylor series 
expansion in the neighborhood of the stagnation point where the second deriva- 
tive is obtained from the radius of curvature found from equation ( 13 ): 


z 



- (1 + 


m, 


S t ' 



(14) 


Equation (l4) is assumed to give a close approximation to the value of z 
over a range of s'/ s £ to a matching point value that remains to be 
determined and depends upon the slope of z at the sonic point. 

Sonic point.- Experimental pressure distribution data for flat-faced, 
round^cornered ^bodies indicate that the slope of the pressure distribution 
curve was inversely proportional to the radius of the body at the point where 
the pressure attained sonic value. The slope or derivative of the pressure 
distribution curve at the sonic point can be estimated quantitatively on the 
basis of an approximate theory discussed in appendix D. The result in natural 
coordinates is 


d 



d 






( 15 ) 


The slope of z at s'/s^ = 1 (sonic point) is 
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in transformed coordinates. 

Matching of curves .- Equation (l4) is assumed to be valid over a range 
of S V S 4 from zero to s^/s^. At s^/s* the curve defined by equation (l4) 
is tangent to a second curve valid over the range s^/s^ < s'/s] < 1. The 
equation of the second curve is assumed to have the same form as equation (l4): 

s = ^ 1 - tj) - c C 1 - f j) (it) 

The independent variable 1 - (s f /s^) is used for convenience. 

Equations (l4) and ( 17 ) are equated at s^/s^. The derivative of equa- 
tion (l4) with respect to s*/s £ is equated to that of equation ( 17 )* Two 
independent equations are then made available to solve for the constant C in 
equation ( 17 ) and the point of tangency, s^/s^. The results are 


ft _ K - m st 
s * m * + m st - 2A 

where A is the absolute value of the coefficient of (s ! /s^) 2 in 
equation (l4) and 



(l8a) 


( 18 b) 


The curves given by equations (l4) and (17) are plotted over their 
respective ranges of applicability in the Z - (s ! /s^) coordinate system. To 
transform the results into the natural coordinate system, p/p s t versus s/s^, 
it is only necessary to superpose the natural coordinates onto the 
z - (s'/s 1 ) curves. This is automatically accomplished with double coordi- 
nate paper (fig. 10). 

The above procedure cannot predict pressure distributions if inflections 
in pressure occur. Such conditions can occur for a large -angle cone with the 
apex blunted by a small radius, R^. However, this case is excluded and the 
method is applicable whenever s £/ s * > 0 (eq. (l8a)). 
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NUMERICAL EXAMPLES 


Shock Solutions 

The following numerical examples illustrate how the shape of the bow 
shock and its proximity to the body are determined for typical vehicles from 
each of the classes considered. 

Example A. Ellipsoid : a/b = 0.25, 7 = 1*4, M = 10 

At M = 10 , the normal shock-density ratio for air is P 2 /p ± = 5*71 or 

Pi/ P2 

sin 0 

*b = 

The values of 4>„ and 0^ are 4^ = sin -1 (0.825 sin 40.4°) = 32*3° and, 
e * = 23.8° (fig. 6 (a)). 

Equation (8) now reduces to 



0 _ A + 16(0. 632) 2 - 1 

x s 

5s = 

16(0.632) 

+ y* 

y* 




L ^ 

R s 


= 0.155 + Z* 
y* Aq Aq 

R s R s 


= 0.175 thus G = 0.116 (fig. 1), A*/y* = O.325 (fig. 2), 

*/sin = 0.825 (fig. 3; a/b = 0.25), 4 = 40.4° (fig. 4), and 

(b/a) 2 = l6 (by definition). 


The value of A q /R s is not known and cannot be determined with equation (lb) 
since the value of R s /Rk is not known. If A q /r s is approximated by 
A 0 /R s fits G = 0.1l6, the equation for R s /y* becomes 


«s 

y* 


01122 + 

0.116 


1 *£. _ 
0.116 y # 


1.336 + 8.621 

y* 


The above equation is now plotted on the coordinates of figure 8 (a). The 
equation is a straight line with the ordinate R s /y„ = 1-336 at x g /y^ = 0 and 
has a slope of 8.621. This line intersects the curve, 0^ = 23.8°, at 

RsA*- 2 -98- 

A trial value for R s /Rb may now be determined from equation 7(a): 
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(a) Ellipsoid. 

Figure 8.- Numerical application of nomograph. 


y. 


* 


0-632 


R b Jl + 16(0. 632) 2 


= 0.232 


and 


R R v 

= — -*• = 2.98 x 0.232 = 0.691 


The value of A 0 /R g ^ from equation (lb), is 


A 0 _n/i + 4 ( 0 . 116 ) ( 0 . 691 ) - 1 


Rc 


2(0.691) 


= 0.108 


which is used to reevaluate R /v : 

s' ■'* 


5s 

y* 


0.155 

0.108 


1 x a , x a 

-4-3 -§■ = 1.435 + 9-259 -§■ 

0.108 y* y* 


and is found to give the value R s /y^ = 3-15 (fig- 8(a)) as the second trial 
solution. New values of R s / R^ and A q /R s are calculated: 

^ = 3.15 X 0.232 = 0.731 
R b 
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and 


Aq = s/1 + 4 (o. 11 6) (0.731) - 1 = o lo8 
R s 2(0.731) 

No change in A /R g has occurred so convergence to a solution has been 
achieved. The value for B s is determined by the intersection point of the 
final equation for R s /y* with the curve, 9 = 23.8° (fig. 8(a)). The final 
results are: 


B s = 4.8 

= 0.108 
R s 

and 

-2 = 0.079 
R b 

Example B. Flat -faced body with rounded corner: r ^/ r rri = 0*25; 7 = 1.4, 

M = 5 

At M = 5, the normal shock -density ratio for air is p 2 /p ± = 5*00 or 
Pl /p 2 = 0.200; thus G = 0.132 (fig. 1), A*/y* = 0-372 (fig. 2), 
sin $ # /sin = 0.712 (fig. 5), and <£*.,_ = 4l.l° (fig. 4). The values e = 0° 

and = 00 are appropriate to a flat -faced body. 

The values of and are §» = sin" 1 ( 0.712 sin 4l.l°) = 27*9° and 
e„ = 21.8° (fig. 6(a)!. 

Equation (10) is now evaluated. (Note that since R^ = oo, equation (l) 
gives A q /R s = G = 0.132.) 

0 372 - o-o + o-o + 0-29 0 (1 - 0 .884) fs 

Rs _ J 1 - 0.230(1 - 0.468) ^ 

y* “ 0.132 0.132 

R« x e 

— = 2.568 + 7.576 — 

y* y* 

The line represented by the above equation has the R g /y^ ordinate value 
of 2.568 at x s /y^ = 0 and a slope of 7*576. This line intersects the curve, 

9 * = 21.8°, in figure 8(b) at R g /y^ = 3*77 and also gives B s 8.0. No 
iteration is necessary since Rg/R-^ = 0 remains fixed. The shock centerline 
radius and standoff distance may be related to the body radius, r , with 
equation 9 (a): 
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y * 


(b) Flat faced body. 

Figure 8.- Continued. 

A = l - 0.250(1 - 0.468) = 0.867 


then 

^■ = — = 3.77 x 0.867 = 3.27 

r m y* r m 

and 

A 0 A n R Q . 

— = — — = 0.132 x 3.27 = 0.432 

r m R s r m 

The final results are B s = 8.0, R s /r m = 3-27, and. A 0 /r m = 0.432. 

Example C. Spherically blunted 65° cone with rounded corner: 

R b/ r m = 0,167 ' r c/ r m = 0,o8 3> 7 = M = 5-2 

At M = 5.2, the normal shock-density ratio of air is p a /p x ~ 5*064 or 
Pi/p 2 = 0.1975; thus G = 0.1305 (fig- 1), A*/y* = 0.368 (fig. 2), 
sin €>*/sin $*•, = 0. 500 (fig. 5), and 3> = 41.05 (fig. 4). The value 

e = 90° - 65° = 25°. 

The angle = sin -1 (0.500 sin 4l.05°) = 19*2°. However, this value 
for is smaller than e = 25°> the angle of tangency between the corner 

radius and the conical surface. Therefore, the sonic point is taken at the 
point of tangency: 
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and has a slope of 7 . 663 . The intersection of this line with the curve 
0 * = 21.4° gives R s /y* = 1*3^- 

A trial value for Rg/R^ may now be determined. Equation (9a) gives 

A = 1 _ 0.083(1 - 0.423) = 0.9519 

r m 

and 

Rs 
R b 

The value of 
for R s /R b , is 


= 4 i + 4(0.130$) (7.651) - 1 = 0 o807 

R s 2(7-651) 

Equation (10), reevaluated with the new value for A 0 /R s , is 


= A Z* = 1.34 0,9 ^ 1 2 = 1.34 x 5*710 = 7*651 
y* r mRb °- l6 T 


A 0 /R s ^ determined with equation (lb) using the above value 


Rs = -0.080 y» 
y* 0.0807 0.0807 


-0.991 + 12.391 


1 

yl 


The intersection of the above line with the curve = 21.4° in 
figure 8 (c) gives R s /y^ = 1*83* The adjusted values for R g /R b and A 0 /R s 
are 


then 


R 

-1 = 1.83 X 5.710 = 10.449 

R b 


Aq = s/l + 4 (0.1305) (10.449) - 1 = 
Rg 2(10.449) 


Rs = -0.080 + y* 
y* 0.0737 0.0737 


-1.085 + 13.569 ^ 


The intersection of the above adjusted line with the 0^ = 21.4° curve 
now gives R s /y^ = 1-95* Another iteration is performed for Rg/R^ and A q /R s : 

5s = 1.95 X 5-710 = 11.134 
R b 
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and 


A o _ n/i + 4(0.1305) (11. 134 ) -l_o 0700 

17“ 2(11.134) f 


Equation (10) is again adjusted to give 


s _ 


-0.080 


2s 

y* 




y* 0.0722 0.0722 


= -1.108 + 13.850 -2. 




The value Rg/y* = 1*97 is now determined from figure 8 (c). No large 
change in R s /y* has occurred, indicating that convergence to a solution has 
been essentially accomplished. The value B s = -2.5 is also determined from 
figure 8 (c). 


The final results are: 


B„ = -2.5 


R 

s. = 1.97 x 5-710 = 11.249 
% 

_ 1 - 4(0.1305)(11.24 9) - 1 


Rc 


2(11.249) 


= 0.0721 


= ^s 2* = 0.0721 X 1.97 x 0.9519 = 0.135 




R s y* ■‘•m 


Example D. Spherical -faced body with rounded corner: Rb/ r m = 2.4, 

r c/ r m = 0.10, equilibrium flow with p 2 /p 1 = 20 

The air is assumed to behave as a perfect gas, and since the Mach number 
is large, the asymptotic value, 7 = [(pg/p^ + l]/[(p 2 /p 1 ) - 1] = 2l/l9 = 1.105* 
characterizes the thermodynamic properties of the air. 

Since no body shock parameter charts are given here for gases with 7 
other than 1.0 and 1.4, interpolations are necessary. 


(The effect of 


The value G = 0.037 I s found from figure 1 at p 1 /p 2 = l/20 = 0.05* 

y on G is negligible at small values of P ± /p 2 " ) The values 
A*/y* = 0.117 is found from figure 2 by linear interpolation between the 
curves for 7 = 1.0 and 7 = 1.4 at p x /p 2 = 0.05* The value 
sin $*/sin = 0.53 is found from figure 5 at r c /r m = 0.10. The value for 
is found from figure 4. Note that there is no great difference between 
the values of <5*! for gases with 7 of 1.0 or 1.4 at the same shock-density 
ratio. The difference in between gases with 7 = 1.1 and 1.0 should be 

negligible, particularly at small values of p 1 /p 2 . Accordingly, the value 
$**1 = 36.8° found from figure 4 at p 1 /p 2 = 0.05 and 7 = 1.0 is assumed valid 
for 7 = 1.105. 
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The value of e for a spherically blunted vehicle of radius R^, tangent 
to a round corner of radius r c , is given by 


r m ( 1 _ *A 

R b \ r m/ 

= sin" 1 

274 (1 - 0.10) 



1 _ 


, 0.10 

L r m R b _ 


L 2.4 J 


23° 


The value 0. 


larger than 
€ = = 23 c 


* 




= sin“ 1 (0.53 sin 36*8°) = 18. 4°. Since the angle e is 
is taken as the sonic point inclination angle; therefore, 


The value 0^ - 0_* o = 4.2° is found from figure 6(b) (p 2 /p x = 20 and 
= 23 0 ). The curves of figure 6(b) are independent of 7. The value of 
0* o , given by equation (c4), is 


0 = tan -1 — = tan" 1 ^*^7 = 17*5° 

*° A* 0.117 


The value of 0^ is found by addition, (0^ - 0^. o ) + 0*o = 17*5° + 4.2° = 21.7°. 
Equation (10) is now evaluated: 


y* 


0.117 - 


0,424 + 2,4 ‘ 0 > a92Q 7 + °- 10 (~ o.92o 91 - 0,920 


1 - 0.10(1 - 0.391) 


Rc 


y* 

% 

r<= 


or 


R s -o.o84 x s /y* 


y* A 0 /R s A q /R s 

If A o/ R s is approximated by G = O.O37, the equation for Rg/y^ becomes 


x c 


+ -J * — = -2.270 + 27.03 — 
y„ 0.037 0.037 y* 


In this case, as in example C, the line representing equation (10) 
crosses the x s /y^ axis at Xg/y^ = 0.084, as is apparent above. The slope 
of the line is 27. 03. The intersection of the line with the 0* = 21.7° curve 
in figure 8(d) gives the value R s /y^ = 2.82. 
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2k 

y# 


(d) Spherical-faced body. 
Figure 8.- Concluded. 


Since the sonic point is at the tangent point of and r c , 
y # = R b sin e , or y^/R b = sin 23° = 0.391* The value for R g /R b is then 


5 §. = !§. Z* = 2.82 x 0.391 = 1.103 

R b y* R h 

The value of A 0 /R q > determined with equation (lb), is 


io = n/i + 1(0.037 ) (l.~103) - 1 
R s 2(1.103) 


0.0356 


Equation (10) is reevaluated: 


h. = - 0.084 + y» 
y* 0.0356 0.0356 


-2.359 + 28.09 


y* 


A new value of R s /y* = 2.93 is found from figure 8(d). Since the value 
of Rg/y^ has not changed greatly, convergence to a solution has been 
essentially achieved. The final values are 


R s ~ 2 

- 2 - = 2.93 x O.391 = 1.146 

R b 
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and 


— = 0.0356 


^2. = fk = 0.0356 x 2.93 x 0.9391 = 0.098 
r m R s y* r m 


Pressure Distribution 

An example calculation for the pressure distribution over the forward 
face of a blunt body is presented in the following paragraphs. Only one exam- 
ple is given since the calculative procedure, with obvious minor variations, 
is applicable to most bodies for which the shock solution has been determined. 

The method is not applicable to certain cases of large -angle cones with 
the apex blunted by a small radius, R^. Criteria for applicability of the 
method are given in the Analysis section. 

The example calculation applies to the round -cornered, flat-faced body at 
Mach number 5 tor which the shock solution was calculated in figure 8 (b). The 
pertinent shock and body parameters are r c / r m = 0-25* = 27- 9 °; 

A o/ r m = pj p 2 = 0 . 20 , and 7 = 1.4. 

First, the numerical values for the coefficients of s T /s_£ and (s^s ^) 2 
of equation (l4) are determined. The values for m^, s */a q , Aq/R-j^, and v 
are required. The quantity m T ^ represents the slope of the line drawn from 
p/p s t = 1 to p^/p s ^ as indicated in sketch (f). In the present example, the 
gas is air for which p/p s t ~ 0 . 528 ; and 


m» = 1 - — =1 - 0.528 = O.U72 
st P s t 

The length, s*, extends from the stagnation point at the vehicle center- 
line to the sonic point at location 0 ^. on the rounded corner of radius r c . 
For the present example, s^ has the value 


s * = “ r c) + (^30) 

which, expressed in terms of A Q , is 

£c\ + Ic 
V W 57-3 r m A o 


s * _ r m 


A o 


or 
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The value for Aq/R^ Is 0 since = oo for a flat face. The value 
v = 0.150 is found from figure 9 at p 1 /p 2 = 0.20 and 7 = 1.4. 



0 .10 .20 .30 .40 .50 

P&2 


Figure 9 * - The v function. 


or 


Equation (l4), written with numerical coefficients, is 


z = 0.472 - (1 + 0.472 2 ) 3/2 [0.150 X 2.018(1 + 0.0) ] 2 

s * 
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Next, the numerical values of the coefficients of 1 - s, / s * an & 

[1 - (s l /s^)] 2 of equation (17) are determined. The values for s^/r c , m£, 
s 4/ s *> anc ^ C are required. 

— = 1“ f* = _L_ 0.432 X 2.018 = 3.487 
r c r c r m A 0 0.25 

Equation (l6) is then evaluated for 


mi 


tan 




tan(cot"‘ :L 0.472 - cot" 1 3-487) 


= tan (64.7° - 16.0°) - 1.144 


The value of A required to calculate the value of the coefficient C 
has already been determined and is the coefficient of (s ! /s_£) 2 in the equa- 
tion for z, that is, A ^ 0.126. The value of the tangent point location 

s’/s’ is also required to evaluate C. Its value (given by eq. (l8a)) is 
p -x- 


ft _ 1.144 - 0.472 

~ 1.144 +0.472 - 2 x 0.126 


0.4927 


The coefficient C (given by eq. (l8b) ) may now be evaluated: 


0.472 + 1.144 - 2 x 0.126 x 0.4927 1 , 71 

2(1 - 0.4927) 


Equation (17), written with numerical coefficients, is 


z 


'■ m G - ij) 


1.471 1 



Values for z 
equation (l4) are 


in the range 0 < s T /s_£ < 0.4927 computed with 


z 


Mlllllllllllllllllll 


s, Ml 

0.0 0.0 

o.i o.o 459 

0.2 0.0894 

0.3 0.1303 

0.4 0.1686 

0.4927 0.2019 

Values for z in the range 0.4927 < s '/ s * <1-0 computed with 
equation ( 17 ) are 

s*/ s 4 z 


0.4927 

0.2019 

0.6 

0.2222 

0*7 

0.2108 

0.8 

0.1697 

0.9 

0.0996 

1.0 

0.0 


The above sets of values are shown plotted on the Z - (s'/s^J) coordi- 
nates of figure 10. The resulting curve represents the pressure distribution 
on the superposed (p/p st ) - ( s / s *) coordinates. 



Figure 10.- The Z - (s'/ s i) coordinate system with an example pressure-distribution solution. 
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COMPARISON OF EXPERIMENTAL AND PREDICTED RESULTS 


The validity of the method was assessed on the basis of the experimental 
results of references 3, 4, 6 , 7 , and 8. These data are not extensive but do 
include representative cases among the classes of body shapes considered in 
the previous sections . 


Shock Shape 

A comparison of predicted and experimental shock traces for various types 
of blunt bodies is shown in figure 11. 



Figure 11.- Comparison of experimental and predicted shock shapes. 
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The geometric parameters of each body are indicated along with the blunt - 
ness, B s , of the associated shock and the Mach number of the free stream. 

Good agreement between experimental and predicted shock shapes is shown. 

In all cases, the shock traces appear circular. However, the shock 
bluntness, B s , has a significant effect on the centerline shock standoff dis- 
tance. To illustrate this effect, an element of a circular shock is shown on. 
the body axis of symmetry for the cases, r c /r m =0.5 and M = 10.53* and ^ or 
the conical body (example C), M = 5*2. The circular shock satisfies the value, 
A */y*, appropriate to the free-stream flow and the value, 0*, appropriate to 
the body geometry. An appreciable error is evident in the shock standoff dis- 
tance at the centerline for the circular shock. This error would increase at 
larger values for the normal shock-density ratio. 

1.0 


.8 

P/Pst 

.6 


.4 





Figure 12.- Comparison of experimental and predicted pressure distributions. 
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The good agreement between predicted and experimental shock shapes veri- 
fies the correlation curve of A*/y* (fig- 2) and the value, 0*, determined 
from figure 5 - The good agreement also confirms the relationship between Q* 
and 0^ (fig- 6) since the value of B s , which strongly influences the center- 
line shock standoff distance, is sensitive to the correct value for 0_*. 


Pressure Distribution 

Predicted and experimental pressure distributions over the forward face 
of various blunt bodies are compared in figure 12 - The results of the present 
method and the theoretical results of reference 1 for the sphere at M = oo 
are also compared. 

Agreement between the predicted and experimental results is usually 
within the scatter of the experimental values. The predicted values are 
somewhat higher than the experimental values for the sharp-cornered, flat- 
faced cylinder at M = 3 - 55 - Values predicted by the method of reference 3 ; 
which is applicable for the two sharp-cornered vehicles shown, are also indi- 
cated in figure 12; those predicted by the present method are in better accord 
with experiment. 


CONCLUDING REMARKS 


A method was developed for predicting shock envelopes and pressure dis- 
tributions for a variety of blunt bodies at zero angle of attack. The method 
is restricted to those cases in which the bow shock is detached from the body 
and the flow over the forward face of the body is subsonic. 

The method is based on correlation functions which relate the shock 
standoff distances at the stagnation and sonic points to the body geometry. 
These correlation functions were developed from the perfect gas solutions of 
reference 1 and depend primarily on the normal-shock density ratio modified 
to a small degree by the specific heat ratio of the gas. Since the effect of 
the specific-heat ratio is small, the present method should give adequate 
solutions for the equilibrium flows of real gases. 

Predicted shock envelopes and pressure distributions were compared with 
experimental values for air flows in the Mach number range of 3-55 to 10.53 
for a variety of body shapes. Satisfactory agreement between predicted and 
experimental values was found for both shock shapes and pressure distributions. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035 , Nov. 13, 1967 
129-01-08-01-00-21 
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APPENDIX A 


SONIC LOCATION ON ELLIPSOIDS 


The sonic location on ellipsoids is found by the Busemann solution for 
infinite shock-density ratio. 

The shock-body geometry in the vicin- 
ity of the sonic point of an ellipsoid is 
shown in sketch (g). The shock-layer 
thickness, cr, is measured normal to the 
body or in the direction inclined at the 
angle,, with respect to the free-stream 

direction. The curve of radius, ft, repre- 
sents the streamline passing through the 
shock at a distance y from the axis of 
symmetry, at which point the shock incli- 
nation Is 6. The Busemann assumptions 
are made that the momentum tangent to the 
shock is preserved along streamlines 
within the shock layer, and that the veloc- 
ity also remains constant along a stream- 
line as indicated. In the limit of a van- 
ishingly thin shock layer, the velocity 
and the stream radius, ft, becomes equal to the local 



Sketch (g) 


is V sin 0 = Vi sin <& 
body radius, R^, at the sonic point. 

A mass -flow continuity equation involving the differential annular flow 
tubes at the shock and within the shock layer, relating dr] and y, is written 


2irp 1 V 1 y dy = sin 0 dr] 


or 

p-,y dy 

= dri 

P 2 y* sin 0 


The values of y and <3> for an ellipsoid are related by 


(Al) 


y _ sin $ 

R b n/ ( B b - 1 ) s in 2 i> + 1 

Equation (A2) is differentiated with respect to $ with the result 

dy _ co s $ d$ 

R b [(B^ - l)sin 2 $ + i] 3/2 


(A2) 


(A3) 


3 ^ 



The variable y is eliminated from equation (Al) with equations (A2) and 

(A3) : 


r -,2 

Pi R b 


cos 4> dO 


[ (B-J-, - l)sin 2 0 + 1]* 


= drj 


(A4) 


The differential pressure, Ap, across the shock layer at the sonic point 
(0 = 0^) is given by 

P^d-p r a p 2 Vf sin 2 0 dr] Pi^iRb P^* sin 2 0 cos 0 do 

AP ^ ^ o R b* y* R b* o [(Bt, - 1) sin 2 0 + l] 2 


The radius, Rb*-> of an ellipse at location is related to the centerline 

radius, R^, by the identity 


Rb t(Bb - ^sin 2 + h 


y* R b* 


sin 0 


(A5) 


* 


The above integration was performed and^ with equation (A5 ) > the result 
was arranged to give the following equation for the differential pressure 
across the shock layer at the sonic point: 


Ap [(B b - 1) s in 2 <&* + l] 2 J 

^ 1 sin -1 

v/ B b - 1 sin 0* 

sin <&* j 

P o.Vf 2(B b - 1) sin ^ 

ye b - 1 

J(B b - i) s in 2 €>* + 1 

( B b - 1 1 s i n 2 0* + 1 J 


(a6) 


An independent equation giving this pressure differential for a gas with 
7 = 1 at infinite Mach number is the Newtonian result 


^2 = U “ Sin2 

Pl V l 


(A7) 


' St 


Equations (A6) and (A7) are equated to give 


[ (B b - 1) s in 2 0* + l] 2 J 

f 1 sin-l 

•J B b - 1 sin 0* 

sin 

2(B b - 1) sin 0* 

k/B b - 1 

s/(B b - 1) sin 2 4>* + 1 

(B b - l) sin 2 0* + 1 J 


(1 - sin 2 <!>*) - (2) 
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TRE 0*. - FUNCTION 


The relationship between 0* and was derived from the results from 
reference 1 for a sphere, along with the results found from the shock solu- 
tions of a series of increasingly (elliptically) blunted bodies. The solution 
for elliptically blunted bodies is developed in this section. 




Sketch (h) 


Sketch (h) depicts a body that is elliptically 
blunted, up to at least the sonic point at y # 

and its associated shock of bluntness B g . 

By inspection. 


A = A. + 
* o 


The values of x s and x-^ are related to y 
through the bluntness parameters, B s and Bt,, 
follows : 


Xs B c 


1 - /I - Be 


J b Vr 


The values of x s and x^ given by equations (B2) and (B3) are substituted 
into equation (Bl) and the result divided by R g to give 


+ R b 

y* R s Rs ®b^'s 


1 - /I - B> 


s z* 

V*b Rs, 


-L 1 - 1 - B S 

B s J S W 


Equation (4) is then solved for y*/R s - The angles 0* and 4>* are then found 
from the following equations, which are derivatives of equations (B2) and (B3), 


tan 0^ = 


11 - B * & 
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tan = 


A 


y* 

R s R’d 



( 6 ) 
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APPENDIX C 


SHOCK BLUNTNESS B g AND ANGLE 


0* o FOR A DISK 


On a flat disk the centerline body radius R^ 


oo and equation (l) gives 


Ao 

R s 


= G 


(Cl) 


The equation of the shock is written in the form 




(C2a) 


where x g is the streamwise distance from the 
shock apex to the location on the shock at the dis- 
tance Av in the free-stream direction from the 
corner of the disk, as indicated in sketch (i). 

The angle 0* o is the shock inclination; the 
length x s has the value A Q - A*. Substituting 
this value for x s in equation (C2a) and with 
R s = A q /G gives 


B 


s 



1 _ 1 ^ 

A* G tan 0 *q/ 


(C2b) 


which is written in the form 


(i _ B ^ = At L_ 

y* \ G s / y* \ A * tan 0 . 



(C2c) 


A solution to equation (C2c) is given if 


and 



tan 0* o 


_G_ 

Ao 

y* 


(C3) 

(C4) 
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APPENDIX D 


PRESSURE DISTRIBUTION 


The pressure distribution at the stagnation and sonic point locations of 
a body are developed as follows. 


Stagnation Point 

In reference 2 it is shown that the stagnation point velocity gradient is 


d V = Vi / AA f 

ds A 0 V R b/ / 2 


(Dla) 


The function, f, of reference 2 is related to the function G as follows: 


Thus, 



dV = Vi 1 + \ £i_ L _ (Pz 

ds A q \ / ^st - \Pi 


1 G 


(Dlb) 


(Die) 


If the velocity, V, is assumed to be linear with s in the vicinity of the 
stagnation point, equation (Die) may be integrated and normalized with respect 
to sonic ( )* values to yield 


V 


Vi f* 1 + V) Pi_ 

V* Ad \ ^b/ p st 


1 



(D2) 


The relationship between pressure and velocity of a perfect gas, written 
in a form valid for small velocities, is 


Z 


Pl . 7 ■ 

1 

(VyP 

1 7-1 

« 1 _ 7 

AT 

L 7 + 

7 

W J 


7+1 

\vj 


Equations (D2) and (10) are combined to give 



( 10 ) 


( 13 ) 
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where 


v 


7 V l 

7 + 1 V* P st 


1 




G 


The function v was determined for gases with y - 1.0 and 1.4 as a function 
of the normal shock-density ratio, P 1 /p 2 (fig- 10). 


Sonic Point 

Newtonian theory gives for pressure distribution. 


P 


P 


st 


cos 2 0 


The derivative of the pressure distribution with respect to $ is 


d 



= -2 sin $ cos $ 


(d4) 


(D5) 


The pressure gradient at the sonic point on a round corner of radius, r c , is 
found by eliminating <S> between equations (d 4) and (D5): 



or, since, generally, 0.5 < P /p a+ . < 0.6, 



Although equation (15) was based on approximate theory, it generally agreed 
well with available experimental data. 
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